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\ This paper focuses on the isospin properties of the asymmetric nuclear-matter 

^ ■ hquid-gas phase transition analyzed in a mean-field approach, using Skyrme effec- 
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tive interactions. We compare two different mechanisms of phase separation for 
low-density matter: equilibrium and spinodal decomposition. The isospin properties 
of the phases are deduced from the free-energy curvature, which contains informa- 
tion both on the average isospin content and on the system fiuctuations. Some 
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^ . implications on experimentally accessible isospin observables are presented. 
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I. INTRODUCTION 

O I The nuclear interaction is a mean-range attractive, short-range repulsive force. Because 

C ■ of this, nuclear matter is expected to present a phase transition of the liquid-gas type 

m. 



The nuclear-matter phase diagram has been studied in detail in the case of symmetric mat- 



^ I ter, i. e. a system with the same number of protons and neutrons [3| . The case of asymmetric 



a. 

matter involves an additional degree of freedom, isospin. As a result, its study is more com- 
plex and allows to observe specific effects such as the phenomenon of isospin fractionation 
^. The knowledge of these properties is of great interest for the physics of neutron-rich 
systems. Such systems have become the subject of many investigations in the recent nuclear- 
physics literature. Exotic-nucleus transport properties and multifragmentation are expected 



to provide constraints on the nu_c 
dence of the symmetry energy 



ear equation of state, in particular on the density depen- 



[5, y, Q]. In astrophysics, both neutron-star structure and 
supernova dynamics are influenced by phase transitions of neutron-rich matter in a large 
interval of temperatures and densities ja, 0]. All these different phenomena involve excited 
matter at baryon densities lower than normal nuclear matter density : this corresponds in 
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the phase diagram to a region of instabihty with respect to phase separation. In the case of 
neutron-star crusts the system hfetime is such that one can safely modehze phase separation 
in terms of equihbrium. However this is not the case for exploding supernova cores, and even 
less for heavy-ion collisions. For such systems with a finite lifetime, the out-of-equilibrium 
)henomenon of spinodal decomposition may be the leading mechanism for phase separation 



1Q|. 

For the work presented in this paper, we have studied nuclear matter in a mean-field 
approach, using Skyrme effective interactions. Three different parameterizations are used. 
The more recent one, Sly230a is used as a reference since it has been fitted to exotic 
nuclei and neutron-matter properties, which makes it particularly well-adapted to the de- 
scription of neutron-rich matter. Two others are used as a comparison : SGII [12], which 
has been shown to reproduce isospin effects studied through giant dipole resonances isl . 
and SIII [l^, which is one of the first Skyrme parametrizations and is not constrained for 
asymmetric matter. Contrary to Sly230a and SGII, SIII is certainly not a realistic param- 
eterization; we have included it in our comparison since it is still used in some transport 
codes because of its simplicity [5|]. 

In section II, we present the phase diagram of nuclear matter, which is a representation 
of phase separation at equilibrium. Then the phase properties are considered, comparing 
the case of equilibrium and spinodal decomposition. In section III, we study the direction 
of phase separation, leading to the phenomenon of isospin fractionation. The consequences 
of our findings on the possible extraction of the symmetry-energy coefficient from isoscaling 
data 0,9 are explored in section IV. In section V we discuss fluctuations in the system 
isotopic composition. 



II. PHASE DIAGRAM 



Gibbs statistical equilibrium is defined by the maximization of the system Shannon en- 
tropy S = —Tr(DlogD), where D is the many-body density matrix of the system [l5|. Let us 
consider S represented for a homogeneous system in the space of observables. If S presents 
a convexity in this representation, it can be maximized by phase mixing, which corresponds 
to a linear interpolation between two points of the space of observables. The coexistence 
region is then defined as the ensemble of points for which entropy can be enhanced by such 
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an interpolation. If the system is located at a given point inside coexistence region, there 
is a unique decomposition in a couple of two other points that gives the best maximization 
of entropy : they define the two phases at equilibrium. This couple is determined by Gibbs 
construction [16|, which consists in building the concave envelope of the entropy surface by 
linear interpolation in the space of observables : phases at equilibrium then belong to the 
same tangential plane [o. \v\. 

This geometrical condition results in the well-known equality of all intensive parameters. 
Let us consider a system occupying a volume V described by a set of observable densities 
{ofc} = {{^k)/V} controlled by the associated intensive parameters {A^}, related to the 
entropy density s = S/V hj : 

Xi = daisiiak}) (1) 
The pressure P is related to the entropy by : 

PP = s-J2 ^kak (2) 

k 

Let us now consider two points of the space of observables {ak}^^^ and {afc}*^^-*. Entropy has 
the same tangent plane at these two points if {XkY^^ = {Afc}*-^-* (equality of plane slopes) 
and P'-^^ = P*^^-* (equality of plane origins) : these are Gibbs equilibrium conditions on the 
intensive parameters. 

In the case of asymmetric nuclear matter, the relevant observables are neutron, pro- 
ton and energy densities, {ak} = {p„,Pp,e}; the associated intensive parameters are 
{Xk} = {—l3fJ'n,—PfJ'p,l3} where fig (g = n,p) are the chemical potentials and [3 is the 
inverse temperature. To obtain the phase diagram of nuclear matter, expressions linking 
observables and the associated intensive parameters (i.e. equations of states) are needed. 
We have worked out such relations in the mean- field approach jisl . 

Let us first express the average energy density of homogeneous nuclear matter, e = 
(H) /V. With a Skyrme effective interaction, e is a functional of one-body densities only: 
particle densities pq and kinetic densities Tg = {P'^)g/h^V . We introduce isoscalar and 
isovector densities : 



P = Pn + Pp , r = Tn + Tp 



PZ — Pn — Pp 1 T3 — Tn — Tp 



(3) 



In the case of homogeneous, spin-saturated matter with no Coulomb interaction, the energy 
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density can be written as : 

e = —r + Cop' + DopI + C3p"+' + D^p^pl + Cefjpr + Z^e/zPsra (4) 
where coefficients Cj an Di are linear combinations of the standard Skyrme parameters 

Co = 3to/8 

Do = -to(2xo + l)/8 

Cs = ts/16 

D3 = ~h{2x3 + l)/48 

Ceff = [3ti+t2(4X2 + 5)]/16 

De// = [t2(2X2 + l)-tl(2Xi + l)]/16 

Mean-field approaches describe nucleons as independent particles. Each particle of type 
q is submitted to the mean field Wg defined by the relation 6{H) = Tr{Wg6pg), where 
6pq represents any variation of the one-body density matrix associated with the particle q, 
leading to the expression : 

The mean field is composed of a potential term Ug = dp^e, and a kinetic term involving 
an effective mass m* defined by h'^/2m* = dr^e. 

This determines individual energy levels = 1^ + Ug with the usual quantification 
conditions on the momenta provided by the boundary conditions. In the grand-canonical 
approach, for a given temperature T = these levels are occupied according to Fermi- 
Dirac statistics : 

ng[ji) = 1/ [1 + exp(/3(pV2m; + f/, - /i,)] (6) 

where Pg is the chemical potential of the particles of type q . 

All one-body densities are then given by Fermi integrals. In particular the independent- 
particle partition sum can be expressed in the same way as for an ideal gas: 

^ = lY^l ln(l+e M^^P = 1:^/5-. (7) 

where p'^ = pg — Ug. At this point it is important to recall that the independent particle 
partition sum Zq should not be confused with the mean-field grand-canonical partition sum 
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which is related to the mean-field pressure log Z = (3PV and can be computed from the 
mean-field entropy using \ogZ = S — (3{< H > —fip < Np > —fin < Nn >)• 

For fixed values of p„ and pp at a given temperature, chemical potentials are obtained 
by iteratively solving the self-consistent relations between pq, pq and m*. We can then 
compute the kinetic density and calculate the entropy as a function of densities. Indeed, the 
mean-field solution comes from a variational approximation to the entropy so the mean-field 
entropy Sq provides the best independent-particle approximation to the exact entropy. The 
mean-field entropy can easily be derived from the independent-particle partition sum and 



the average of the mean-field potential {W)o so we get [191] : 



S/V - So/V = \nZo/V + (3i{W)o/V - p^Pn - fipPp) (8) 

which gives for a fixed (3 the entropy s = S/V as a function of neutron and proton densities. 

As we have already mentioned, the construction of the concave envelope of s{e,pn,Pp) 
imphes Gibbs conditions for each couple of phases at equilibrium : f3^^^ = j3^'^\ pn^ = pn\ 
Pp^ = p'p^ and P^^^ = P^'^\ Let us now introduce the constrained entropy defined by 

Sc = -f3f = s- I3e. (9) 

which corresponds to the opposite of the free energy /, divided by the temperature. Gibbs 
construction can then be performed by constructing the concave envelope of Sc in [p^, pp) 
plane at a fixed temperature T = Indeed, /? is imposed to be constant, and the defini- 
tion of Sc ensures that its tangent-plane slopes and altitude give respectively the chemical 
potentials and the pressure : 

{dp,Sc{pn, Pp))p = dp^s{e, pn, pp) = -(3pq (10) 

Sc + PiPnPn + PpPp) = S- (3{e- PnPn - PpPp) = PP (H) 

The introduction of the constrained entropy permits then to reduce the dimensionality 
of the problem since we can work at constant temperature. The constrained entropy is 
represented in figlH at a temperature of 10 MeV. In this figure, a slope has been subtracted 
in order to emphasize curvature properties: the surface represented is s'^ = Sc + PpsP, where 
Ps is the chemical potential of symmetric matter at the transition. The constrained entropy 
is symmetric with respect to the axis p„ = pp, refiecting the invariance of nuclear interaction 
with respect to neutron-proton exchange. 
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T=10 MeV 




FIG. 1: Constrained entropy at T = WMeV, corrected by PfJ,sP in order to emphasize curvature 
properties. It is shown on the left in perspective, and on the right as a contour plot. On the right, 
coexistence and spinodal boundaries are added as thick hues. Outer curve: coexistence. Straight 
Unes Unk selected couples of phases at equilibrium. Inner curve: spinodal. Dots: critical points. 

We observe that the surface obtained presents a region of positive convexity : this is the 
case if we fix any temperature under the critical temperature of symmetric matter, which is, 
in the Sly230a parameterization, Tc = 14.54MeV (for comparison, we get Tc = lAAQMeV 
with SGII and 17.96Mel^ with SIII). The concave-envelope construction performed on this 
surface links couples of phases at equilibrium. The ensemble of these points forms the two 
branches of the coexistence curve, one at lower density (gas) and the other at higher density 
(liquid). They join at two critical points, where the difference between the two phases 
disappears. 

The coexistence curve is the boundary of the coexistence region, inside which a system at 
equilibrium is decomposed into two phases situated on each branch. The region presenting 
a positive convexity is the spinodal region. In this region, infinitesimal density fluctuations 
may lead to an increase of entropy. The system is locally unstable. The spinodal region is 
by construction inside coexistence. Coexistence and spinodal boundaries are tangent at the 
critical points. These curves are represented on the right part of figO The phase diagram 
of the three different Skyrme forces is compared for the same temperature T = 10 MeV 
on the left part of figEl Sly230a and SGII give close results, while SIII shows an atypical 
behaviour, especially, as expected, on the asymmetric parts of the diagram. 
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FIG. 2: Coexistence and spinodal curve in the density plane at T = lOMeV. Left: comparison 
between different Skyrme parameterizations : Sly230a (thick hnes), SGII (thin hnes) and SIII 
(dashed Hnes). Right: illustration of phase separation for a system of proton fraction Z/A = 0.3 
inside instability region (c/ section III). 

III. DIRECTION OF PHASE SEPARATION 



Let us consider a system whose proton fraction Z/A and total density p are such that it 
is situated inside instability region. This is illustrated on the right-hand side of figE] for a 
neutron-rich system {Z/A = 0.3). 

If equilibrium is reached, this system will undergo phase separation according to Gibbs 
construction. We can see that the corresponding phases in coexistence, represented as black 
dots on the coexistence border on the right part of figlSl do not belong to the line of constant 
proton fraction. We can remark that the liquid fraction is closer to symmetric matter than 
the gas fraction: this behavior is a consequence of the symmetry-energy minimization in the 
dense phase. This inequal repartition of isospin between the two phases is the well-known 
phenomenon of isospin fractionation 

In order to quantify the phenomenon of isospin fractionation, we will study the associated 
directions in the density plane. The line of constant proton fraction defines a direction uz/a- 
Let us introduce the isoscalar direction Up, which corresponds to the axis of total density 
P = Pn + Pp- The equilibrium direction Ugg is given by the straight line that joins the couple 
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of phases in coexistence. Isospin fractionation is linked to the difference between u^q and 
uz/A- The equihbrium direction being rotated towards Up, phase separation gives a hquid 
more asymmetric than the gas. 

We now turn to study isospin fractionation if phase separation occurs out of equihbrium. 
SpecificaUy, if the system is brought inside the spinodal region in a reaction too fast for 
global equilibrium to be achieved, its evolution will be driven by local instabilities instead 
of global Gibbs construction. This is spinodal decomposition [lO| : the local properties of 
the constrained entropy curvature determine the development of a spinodal instability into 
phase separation. The information for this study is contained in the free-energy curvature 
matrix given by : 

^ ^ j d'^f/dpl d'^f/dpndpp \ [ dpn/dpn dpp/dpn . ^^2) 
yd^f/dppdpn d'^f/dpl J ydfin/dpp dpp/dpp 
This matrix is defined for each point of the density plane. The lower eigen-value C< 
corresponds to the minimal curvature of / at the considered point. The spinodal border 
is the ensemble of points for which C< = 0. The instability region is defined by C< < : 
in such points, there are directions of abnormal curvature, corresponding to a concave free 
energy. The eigen-vetctor u< associated with the lower eigen-value gives the direction of 
most negative curvature: this is the instability direction. It is represented by the double 
arrow on the example of figj2l where we can see that spinodal decomposition leads to a more 
pronounced fractionation than equilibrium, the denser phase getting closer to the symmetric 
matter. It is important to remark at this point that the instability direction does not give 
a direct quantitative measurement of the final fractionation after phase separation. For this 
to be true, the effect of the evolution of density fluctuations on the isospin composition has 



to be completely neglected. If this is certainly an approximation 27|, it has been however 



shown in explicit transport calculations 10| that the initial instability largely dominates the 
subsequent dynamics. The interest of employing this simplifled approach as complementary 
to more sophisticated transport codes is that, within this approximation, the spinodal mech- 
anism can be directly compared to phase separation at equilibrium. Indeed the equilibrium 
direction can also be expressed as the eigen- vector of a curvature matrix, like the instability 
direction. Let us deflne /^'^ as the free energy at equilibrium, i.e. after correction of the 
curvature anomaly by phase mixing according to Gibbs construction. This construction 
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imposes a positive curvature at each point, and a zero curvature between two phases in 
coexistence. As a result, in the coexistence region the lower eigen- value of the f^"^ curvature 
matrix is zero, and the associated eigen-vector is what we have defined as the equilibrium 
direction, Ueq. In the following, we will denote /° and /'^'^ the free energy before and af- 
ter Gibbs construction, their respective curvature eigen- modes being (C° , u° , C° , u^) and 
{Ct? , , Cy , Uy) . A qualitative impression of how equilibrium and instability directions 
compare through the density plane can be gained from figl3] for the three Skyrme forces. 
In order to follow the eigen-direction of the curvature matrix, we draw curves which are 
everywhere tangential to the curvature eigenvectors. Both directions of phase separation at 
equilibrium (Gibbs construction) and out of equilibrium (spinodal decomposition) remain 
quite close, and essentially isoscalar ^, A comparable degree of isospin fractionation, 
with formation of a more symmetric liquid, is conserved all through the instability region. 
The dominant feature is that the instability direction leads to more fractionation than equi- 
librium. 
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FIG. 3: Global representation of phase-separation directions through the density plane, for three 
different Skyrme parameterizations. Coexistence and spinodal boundaries are represented as in 
figl2j A few selected equilibrium directions are given by straight lines joining the high and low- 
density branches of the coexistence curve (these straight lines are everywhere tangent to the equi- 
librium direction u^^). Phase-separation directions in spinodal decomposition are represented by 
a selection of lines which are everywhere tangent to the instability direction u^. They are drawn 
as full lines inside the spinodal region where they correspond to an instability direction, and pro- 
longated in dashed lines ouside the spinodal, where there is no more local instability. 
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In order to have a quantitative comparison of the different fractionations, it is convenient 
to express the directions of phase separation as slopes with respect to the isoscalar direction 
5p^/5p. This quantity gives a measure of the degree of fractionation associated with phase 
separation. It can in principle vary from zero (in the case of a purely isoscalar phase 
separation), to infinity (in the case of a purely isovector instability, corresponding to a 
separation of protons from neutrons into two phases at the same baryon density). In practice 
6p-i/6p is always a small number [5|,[2^. This is easy to understand considering that i) proton 
and neutron kinetic energies favour equal-size Fermi spheres and ii) the nuclear interaction is 
more attractive in the proton-neutron than the proton-proton and neutron-neutron channels 
in any realistic nuclear physics model (T = contribution), which favours the existence of 
a mixed isospin bound phase. 

The evolution of 5p^/5p (degree of fractionation) with ps/p (system asymmetry) is given 
by figJUfor three different total densities : below, at, and above the respective critical density 
Pc {i.e the total density at critical points). In this representation, the direction of constant 
proton fraction is a diagonal straight line while the isoscalar direction is an horizontal line. 
On this figure, equilibrium and instability direction are compared for the different Skyrme 
parameterizations. 

Phase-separation evolutions verify some expected constraints. First, the isospin invari- 
ance of Skyrme interactions ensures that the graphs are symmetric with respect to p3 = 0. 
It also imposes that all directions become purely isoscalar in the case of symmetric matter, 
i.e. all curves must join at the origin {5pz/ 5p, p^/ p) = (0,0). Another constraint is given by 
the ending points of the curves. At these borders (Ipal = p), the vanishing of one species 
means that its chemical potential at finite temperature goes to — oo, implying an infinite 
curvature of the free energy in the corresponding direction. The minimal curvature u° is 
then oriented along the axis of the remaining species, which is also the direction of constant 
Z/A. 

The first observation is that spinodal and phase-equilibrium directions are very close and 
lie between the isoscalar direction and the constant-proton-fraction one. This means that 
they lead to a rather similar fractionation. In order to get a deeper insight, let us consider 
the various densities in more details. 

Let us consider the case of the axis at p = pc, for which the instability region ends 
at the critical points. Between these extremities, the instability direction induces more 
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FIG. 4: Phase-separation directions 5p^/5p at equilibrium (dash) and in spinodal decomposition 
(soHd) as a function of the asymmetry ps/p for three different Skyrme parameterizations at different 
values of p : below, at and above the associated critical density. The instability direction is 
continued by the lowest curvature direction (dotted line). The dot-dash straight line gives the 
direction of constant Z/A while the double-dot-dash horizontal line gives the isoscalar direction 
^P3 = 0. Small dots give spinodal boundaries and large ones coexistence region boundaries. 



fractionation than equilibrium since it is closer to isoscalar direction. At the critical points, 
where the spinodal and coexistence curves are tangent, both directions of phase separation 
join to become tangential to the coexistence and spinodal borders. This is an additional 
constraint on phase-separation directions which, with the already discussed constraints, 
explains the similarities between the phase-separation directions at and out of equilibrium. 
This additional constraint does not apply if the total density we consider is not pc- We 
observe that for p < pc, a. small inversion occurs around the extremities of the instability 
region. For p > Pc instead, a phase separation out of equilibrium induces more fractionation 
for all isospin asymmetries. 



12 



These observations stand for the three Skyrme forces that we have considered. The 
quantitative differences, which essentially distinguish SIII from the two more recent param- 
eterizations, do not affect the qualitative features that have been described. This is partly 
due to the above discussed constraints. 



IV. FRACTIONATION AND ISOCALING OBSERVABLES 



In a variety of experimental situations, it has been observed that the isotopic-yield ratio 
for two reactions (denoted (1) and (2)) of different asymmetry approximatly follows an 
exponential form |6|: 

Rtl''^^\N, Z) = r(2)(iv, Z)/Y^^\N, Z) oc exp(a„iV + a^Z). (13) 

This phenomenon is known under the name of isoscaling. In an equilibrium interpretation in 
the grand-canonical ensemble, the isoscaling parameters should be linked to the chemical 
potentials and temperature of the systems, aq = {jif^ — jj,^q^)/T. It has been proposed 



to extract from the measured isosca^ 



through the approximate formula 



a 



ing parameters the symmetry-energy coefficient C, 



sym 



an = AC,yU{Z/{A)]l) - [Z/{A%^)/T. (14) 

This relation, derived for the grand-canonical ensemble in the saddle-point approximation 
[7], holds in the hypothesis of equilibrium at the time of fragment formation. This means 
that to extract Csym{Z) from the measured isoscaling parameters, modifications of a„ due 
to secondary decay have to be accounted for and the proton fraction of the most probable 
isobar for each Z should be known at the time of fragment formation . 

Because of the small variation of with the atomic number, an average estimation of 
Csym. can be obtained using for [Z/ the average proton content of the liquid fraction j^]. 

Since even this quantity is not experimentally known, a similar formula can be found in the 
literature, where, however, the proton fraction of the fra gme nting source [Zo/Ao](i) is used 



instead of [Z/ {A)](^iy. this leads to a quantity we call C^y^^ [2l|, l22|, l24] • This approximation is 
supported by the fact that in the zero-temperature limit the proton fraction of the fragments 
converges to the proton fraction of the source {22 1. 

However, in general this does not need to be the case and the approximated C^y^ may 
differ from the value Csym- Iii particular, both in the equilibrium hypothesis at finite tem- 
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perature and in the case of dynamical fragment formation through spinodal decomposition, 
as we have seen, we can expect the phenomenon of fractionation: fragments should not 
have the same isospin ratio as the fragmenting source. This should affect the value of Csym 
obtained from the isoscaling coefficients. 

Our equilibrium nuclear-matter approach does not allow to give quantitative predictions 
on the actual isospin content of fragments obtained out of the different phase-separation 
mechanisms; we can however make some qualitative order of magnitude estimations of the 
possible deformation that isospin fractionation may induce on symmetry-energy measure- 
ments. 

Dynamical transport models 5|, [Zl] indicate that clusterization takes place at a relatively 
well-defined density inside the spinodal region determined by the collision dynamics and cor- 
responding to the first development of instabilities, for which we take a typical value around 
Po/3. Spinodal decomposition is then a phenomenon fast enough for fragment formation 
dynamics to be dominated by the amplification of the most unstable modes at this density 
10l |. This density then determines the phase-separation direction. 

The left part of figl5] displays, similar to figJH above, the directions of phase equilibrium 
and spinodal decomposition as a function of Z/A at a typical temperature of T = 4MeV. 
The case of two systems with global proton fractions corresponding to ^^^Sn and ^^"^Sn 
is shown by the vertical lines. We can see that fractionation at equilibrium is negligible 
in agreement with ref.j22|, and fragmentation of more exotic nuclei should be studied to 
explore this effect. However even for these stable systems, spinodal decomposition gives an 
important fractionation, and the effect would be amplified using more asymmetric systems. 

Once the direction of phase separation is known, the final composition of the fragments 
depends only on their density at freeze-out, when they become separated. We call this 
density psep- The average Z/A of the liquid fraction is represented in the central part of figl5] 
as a function of the density at separation. This density must be between the initial density 
and the saturation value, but its actual value is not known from fundamental principles. It 
is generally assumed that fragments at separation are at normal density psep ~ Po- Recent 



AMD ca: 
realistic 



culations however seem to indicate that lower values p ~ 0.1 /m ^ may be more 
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The asymmetry coefficient Cgym deduced through eq. flHl) is compared to its approximate 
value Cgy^ obtained with the proton fraction of the fragmenting source in the right part of 
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FIG. 5: Influence of isospin fractionation on Csym measurement at T = AMeV . Left side: phase- 
separation directions. Top-left: separation direction at a typical density of fragment formation 
P = Pini = 0.05/m~'^, as a function of the global proton fraction for equilibrium (dashed line) 
and spinodal instability (full line). Vertical lines indicate the two systems Sl~^^^Sn, 
Bottom-left: representation of the corresponding phase separations in the density plane. Right 
side: consequence on symmetry-energy measurement. Top-right : average isotopic composition 
of the liquid- like fraction as a function of fragment density at separation, for the two systems SI 
and S2. Bottom-right: Symmetry energy coefficient Cgym compared to the approximate expression 
obtained neglecting fractionation C^y^, as a function of separation density, for equilibrium (dashed 
line) and spinodal instability (full line). In the case of equilibrium, the liquid position is fixed in 
the density plane: it is indicated by the empty circles. 
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figEl Little change (few % at maximum) is observed in the case of equihbrium, again consis- 
tently with the predictions of the statistical model j22|. However, the stronger fractionation 
effect in the case of spinodal decomposition sensitively decreases, by up to 50% (or more), 
the approximate asymmetry coefficient. This implies that the low values reported for Cgym 
in ref. 2l|] may be alternatively interpreted as a signal of isospin fractionation in a spinodal 
decomposition. 



V. ISOSPIN FLUCTUATIONS 

In section llllt we have studied the direction of phase separation, corresponding to the di- 
rection of minimal curvature of the free energy u<. This direction can be qualitatively linked 
to the average isospin content of fragments issuing from the phase-separation mechanism. 
As we have seen in section IIVI some quantitative predictions can also be obtained if some 
extra hypotheses are made, namely that fragment-formation dynamics is dominated by an 
amplification of instabilities from well-defined initial density and temperature. Within the 
same line of reasoning, the other eigen-vector u> giving the direction of maximal curvature 
can be linked to fluctuations in the system composition, that we now turn to analyze. 

A physical system characterized by a given value of the isoscalar density p at equilibrium 
presents a distribution of ps given by: 

PA.,.3(P3)oce^(^^^"^(^3)+/3^P3) (15) 

where /is = /in — /ip is the isovector chemical potential, and /i = /i^ + /ip the isoscalar 
chemical potential. The constrained entropy s^'^{p^) is deflned by the Legendre transform 
(c/ eq.dH)): 

5f''^(P3) = 5(e,p,P3)-/5e + /?|p (16) 

In the case of an inflnite system, the distribution ( |T5i) is characterized by a unique value, 
P(p3) oc (5(p3 — P3), given by the maximum pg of the function s(?'^(p3) + P^ps- Indeed the 
correspondence between densities and chemical potentials is one-to-one for inflnite systems 
since fluctuations go to zero at the thermodynamical limit. If conversely we consider systems 
of finite volume V, this distribution can be approximated by a gaussian if the system is big 
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(P3-P3)^ 
^2 



enough: 

Pp,,,M^e -i^3 (17) 

where the distribution variance =< p\ > — < which measures the isospin 

fluctuation, is directly linked to the entropy curvature. Indeed if the entropy has a finite 
curvature in ps direction, a development of s^'^{ps) around pg gives: 



^^•^ (P3) + /?f P3 = [s^^^'ipD + P^P^l + ^g^, lp° (18) 

2 1 /^S^sf'^V' 1 /^^'/^ 



.Am^^On , nP3_^0 



(p3-p0)'92sf''^, 



SO that 



This is related to the symmetry-energy coefficient Cgym defined from the free-energy curva 



ture by d'^f /dpi = ^Csym{p), which is independent of pa in the parabolic approximation [29| 

,2 P 



a 



2(3VCsyUp) 



(20) 

The neglect of higher order terms in eg. (1181) is not fully justified for a genuinely finite 
system, where we expect deviations from the gaussian ansatz (fT7|l : it is however consistent 
with our mean-field approximation. Equations f[T7l) . fl2Ul) show that isospin distributions are 
directly linked to the symmetry properties of the effective interaction. This is why isotopic 
distributions measured in fragmentation reactions have been studied by many authors and 
their widths tentatively connected to the symmetry-energy coefficient Cgym. at subsaturation 



density |22|,|26|. 



Let us now consider a system at equilibrium inside the spinodal region at a point (p„, pp). 
Since this situation is unstable for an homogeneous system, it will separate into a liquid and 
a gas fraction. As we have discussed in section IlIIl this separation follows the direction of 
minimal curvature u< of the free energy. The directions of the two orthogonal eigenvectors 
u< and u> define a set of orthogonal eigen-observables (p<,p>) corresponding to a rotation 
of the set (p„,pp) with an angle a = (u<, u„): 

P< = pnCosa + ppsina 

(21) 

p> = Pnsina — PpCosa 
At each point of the system evolution along the axis p<, the inverse curvature along the 
orthogonal axis measures the fluctuation of p> . Since the direction of phase separation u< 
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is essentially isoscalar (a f5i 7r/4), the observable p< is close to the total system density while 
p> fs ^"^^^ = approximately gives the isospin content of the system. In that sense, the 
fluctuations of p> are related to the isospin fluctuations of the system. The actual isospin 
fluctuations can be calculated projecting the fluctuations of p> on the pa direction. 

If the process occurs at equilibrium, phase separation follows Gibbs construction rules, 
which constrain the two phases on the coexistence curve. Thus, the physical equilibrium 
fluctuations for the system are associated with the curvature of the free-energy surface 
f^'^iPmPp) in the direction pcoL (Pcog), given by the tangent to the coexistence curve at the 
liquid (gas) point. 

Since phase separation is a simple straight line, these final fiuctuations are nothing but 
the extrapolation to the coexistence border of the /^"^ inverse curvature in the u> direction, 
calculated at the initial point {pn, Pp) ■ 
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FIG. 6: Inverse curvature of the equilibrated free energy as a function of the isovector chemical 
potential ^3. Dashed lines: curvature taken in the direction ut5 for points inside the coexistence 
region along the axis p = pc- Thick (thin) lines: curvature taken in the direction UcoL (^cog) 
tangent to the liquid (gas) branch of the coexistence curve. Left part: T = 4MeV. Right part: 
T = lOMeV. Lower part: projection on ps axis of the different curvatures considered. 
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The different curvatures associated with the two phases in coexistence are represented 
in figM as a function of the isovector chemical potential fi^ = fin — fJ'p, in order to have a 
common abscissa for liquid and gas. The projections in the ps direction are given in the 
lower part of the figure. The curvature of f^"^ is almost constant inside the coexistence zone 
(dashed lines). Because of the characteristic shape of the coexistence curve (c/ figEI) this 
curvature is reduced (increased) at the liquid (gas) border. It is represented for both phases 
in full lines. This very general result is only due to the free-energy structure of the two 
phases, and we may expect it to stay qualitatively valid independently of finite-size effects 
and phase-separation mechanism. The extreme values of /is correspond to the critical points, 
which are defined by a zero curvature along the coexistence curve. At these points, both 
curvatures in pcoi and pcoc become zero, fourth order terms have to be considered in the 
expansion f|T8|) . and the gaussian approximation f|T7j) breaks down. Since the coexistence 
curve is closed, there is a value of asymmetry for which the coexistence line is tangent to the 
isoscalar direction, leading to a divergence of the ps-projected curvature on the liquid side. 
However, this still happens close to the critical point, where the gaussian approximation is 
not enough. 

Far from the critical point, these inverse curvatures can be associated with a fluctuation 
through eq. (fT9|) . These fluctuations are largely independent of the system temperature 
inside the coexistence zone. They can be strictly interpreted as isospin fluctuations only for 
symmetric matter (ps = 0), where the pcoL and pcoc directions are purely isovector. For 
symmetric matter, the liquid fluctuation appears much larger than the gas one. Because of 
the actual shape of the coexistence line, the isospin part of the liquid fluctuation is a rapidly 
decreasing function of the global asymmetry, as shown in the lower part of figj6l Conversely, 
the gas branch presents a fluctuation increasing with the asymmetry. 

As far as physical fragmenting systems are concerned, fluctuations are affected by con- 
servation laws which are not accounted for in this grand-canonical approach : therefore, the 
fluctuations represented in the lower part of figj6] can only be qualitatively interpreted as 
widths of fragment isotopic distributions. 

Let us now focus on the difference between equilibrium fluctuations and those expected in 
spinodal decomposition. In both cases we will use the same physical picture so that results 
can be qualitatively compared. 

In the preceeding discussion we have assumed that phase separation is entirely deter- 
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le linear 



27|, and 



mined by equilibrium rules. If, out of equilibrium, spinodal decomposition is the dominant 
mechanism, fluctuations in the two dimensional (p„,pp) plane grow in time from the insta- 
bility region until they cause the decomposition of the system. Then isospin fluctuations 
should be calculated as a dynamical variable given by the projection over the isovector axis 
of the two dimensional fluctuations varying in time. This approach is followed in numerical 
codes sl, as well as in simplified analytical models like in ref. [27I, which follows the 
time evolution of unstable modes under the influence of a stochastic mean field in t 
approximation. Both approaches give similar predictions for the isotopic widths 
suggest narrower distributions than in experimental data and in the case of fragmentation 
in statistical equilibrium 26|. 

In these approaches, the initial condition is given by statistical equilibrium inside spinodal 
region of the observables which are not associated with the instability. The initial fluctuation 
of these observables is neglected. However, the inverse curvature along the axis orthogonal 
to the instability direction is never zero, giving rise to an initial fluctuation of p> . We 
now want to estimate this effect. Consistently with the reasoning of section [IV] above, we 
assume that this fluctuation is mainly determined at the onset of phase separation, which 



is justified in ref. 



28|. 



As we can see from figJTl this inverse curvature is not negligible but rather comparable 
to the case of equilibrium. In this figure, we compare the cases of equilibrium and spinodal 
decomposition for the three different Skyrme parameterizations. We consider systems along 
the axis joining the critical points. The abscissa is the isovector density of the global system, 
Ps = Pn~Pp- For each point, we represent the ps projection of the different inverse curvatures. 
For spinodal decomposition, the relevant free energy is /°, and its inverse curvature is 
calculated in direction. At equilibrium, the concave intruder in is eliminated by 
Gibbs construction, and the inverse curvature of the resulting free energy /^'^ is computed 
in u^' direction for the global system, which can be directly compared to the convexity of 
the unstable homogeneous system. 

The main features obtained do not depend on Skyrme parameterization. We can see that, 
inside spinodal region, the convexification of the free energy does not modify drastically the 
curvature in the isospin direction. If spinodal decomposition is the dominant mechanism of 
phase separation, the isospin fluctuation tends to decrease with increasing global asymme- 
try, as in the equilibrium case. A complete dynamical treatment of the problem is needed 
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FIG. 7: Study of isospin fluctuations, compared at equilibrium and for spinodal decomposition 
for three different Skyrme parameterizations as a function of the total isovector density. The 
represented quantity is the inverse curvature of the free energy taken in the relevant direction and 
projected on /93 axis. Global systems are taken along the axis p = pc- The system undergoing 
spinodal decomposition is represented by thin lines (curvature of taken in u2, direction). The 
curvature in vl^ direction is shown by the dashed lines. Dots give the position of critical points. 
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to have quantitative predictions for the final fragment isotopic widths 
the minimal hypothesis of a completely diabatic dynamics that propagates the initial fiuc- 
tuation over the free energy surface, we can expect the absolute isotopic width of fragments 
issuing from spinodal decomposition to be comparable to the case of coexisting phases at 
equilibrium. The high asymmetry region where this trend is inversed is not physical, since 
the gaussian approximation breaks down there for the liquid fraction in the equilibrium case. 

To get more insight on this issue, we can compare the (initial) fiuctuations expected for the 
spinodal mechanism with the predictions of the equilibrium SMM model in the grancanonical 
approximation. Within this formalism the production probability of an isotope composed 
of Z protons and N neutrons reads 



(22) 



Let us denote ^3 = Vp^ = N ~ Z. The isospin variance of the distribution fl22|) is 
i'^AaY ^3 > ~ < ^3 >^= {'^plY- Using the same saddle point approximation we 
used to get eq. ffT^ . this is related to the free-energy curvature by 

1 fd^F\~^ A 



2/3C^jjym(^) 



(23) 
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where Csymi^) is the symmetry energy of the produced fragment. Denoting p*^' the freeze- 
out fragment density, we have Cgymi^) = Csymi.P^'^)- 

In the minimal hypothesis that the initial fluctuations (cr^'^)^ are not amplified (per unit 
volume) by the successive dynamics of phase separation : 



{<)]iyf = {<tiy^ = y^{<S.^ (24) 

the final variance in the case of spinodal decomposition (ca^)/ reads (c/ eq. (|T9l) and (l20l) ) : 



/ .dN2 _ VfPi _ A Pi 

where pi is the density at the onset of fragment formation and pj'^ the freeze-out fragment 
density in a spinodal decomposition scenario. Equations fl23l) and fl25l) lead to the relation : 

«)J _ Csy^ip^'^) P. ^26) 



{o'^aS ^symiPi) pf 



Using the standard representation [29|] Csym{p) = Csym{po){p/ Po^ and considering pj^ ^ 
p'^'i ^ po, we get 



which can be even greater than 1 for an asy-stiff equation of state. 

Considering that in the dynamics of spinodal decomposition isovector density fluctuations 
will also be (slightly) increased from the initial condition 27|], eq. (!27I) confirms that the 
fluctuations associated with the two scenarios are expected to be comparable. 



VI. CONCLUSION 



In this paper, we have presented isospin properties of the phases formed in nuclear- 
matter liquid-gas transition. The characteristics of phase separation are deduced from the 
free-energy curvature properties, / corresponding to the constrained entropy at a fixed tem- 
perature (3. Under a critical temperature, the free energy of the homogeneous system /° 
presents a region of abnormal curvature, where phase separation is favorable. Phase equilib- 
rium is obtained according to Gibbs construction, determining an equilibrated free energy 
f^'^. Phase properties are then deduced from the free-energy curvature matrix, studying 
/'^'^ for a system at equilibrium and /° for a system undergoing spinodal decomposition. 
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The direction of phase separation is given by the eigen- vector u< hnked to the lower eigen- 
value of the curvature matrix. In both equilibrium and spinodal decomposition, it leads to 
isospin fractionation, with a liquid more symmetric than the global system. In the direc- 
tion orthogonal to phase separation, the curvature is linked to fluctuations affecting phase 
composition. This direction being close to p3 direction, it is strongly linked to isospin fluc- 
tuations. Spinodal decomposition is predicted to give a stronger average fractionation and 
at least comparable or even higher isospin fluctuations with respect to a phase separation 
at equilibrium. Isoscaling observables as well as isotopic distributions are expected to be 
sensitive to such effects. If fragmentation can be associated with a spinodal decomposition, 
we qualitatively expect i) low apparent values of the symmetry energy coefficient extracted 
from isoscaling analyses, and ii) isotopic widths comparable or even larger than in the case 
of the statistical model. Finally we would like to stress that the present work concerns an 
idealized study of bulk matter, which is likely subject to significant modification when the 
nuclear surface is taken into account. The pursuit of this issue would be a very interesting 
topic for further study. 
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